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Introduction; Sodium in the lunar exosphere is re- 
leased from the lunar regolith by several mechanisms. 
These mechanisms include photon stimulated desorp- 
tion (PSD), impact vaporization, electron stimulated 
desorption, and ion sputtering. Usually, PSD domi- 
nates; however, transient events can temporarily en- 
hance other release mechanisms so that they are domi- 
nant. Examples of transient events include meteor 
showers and coronal mass ejections. 

The interaction between sodium and the regolith is 
important in detemuning the density and spatial distri- 
bution of sodium in the lunar exosphere. The tempera- 
ture at which sodium sticks to the surface is one factor. 
In addition, the amount of thermal accommodation 
during the encounter between the sodium atom and the 
surface affects the exospheric distribution. Finally, the 
fraction of particles that are stuck when the surface is 
cold that are rereleased when the surface warms up 
also affects the exospheric density. 

In [1], we showed the “ambient” sodium exosphere 
from Monte Carlo modeling with a fixed source rate 
and fixed surface interaction parameters. We com- 
pared the enhancement when a CME passes the Moon 
to the ambient conditions. Here, we compare model 
results to data in order to determine the source rates 
and surface interaction parameters that provide the best 
fit of the model to the data. 

Data: Observations of the lunar sodium exosphere 
were obtained at the McMath-Pierce solar telescope 
with the echelle spectrograph at high resolution during 
November 1998, and January and February, 1999. The 
slit was positioned perp^idicular to the equatorial east 
and west limbs and the north and south poles, respec- 
tively from the surface to about a half lunar radius 
above the limb. The density at the surface and the scale 
height will be determined at these cardinal points and 
compared with models. These data, which have not 
been published previously, and additional data that are 
published, including a sequence in June 1998 [2], 
taken with the same telescope and setup, will be used 
to constrain models of the lunar exosphere. 

Model: The Monte Carlo model follows 10^-10* 
particles on ballistic trajectories from their points of 
release on the lunar surface until the particles are lost 
from the atmosphere using a fourth order Runge-Kutta 
(RK-4) algorithm [1,3). The results presented here use 
the equation of motion including radiation pressure, 
but neglecting the effects of surface charging. The 
simulation space spans fit>m the surfiice of the Moon to 


the Hill Sphere (35 Rmoot), where gravity from Earth 
begins to dominate the motion of the particles. 

An input flux and spatial distribution is assigned as 
appropriate for the source: solar UV radiation for pho- 
ton-stimulated desorption, solar particle flux for ion 
sputtering, and micrometeoritic or meteoritic flux for 
impact vaporization. At the Moon, solar wind flux 
dies off with solar zenitii angle due to the curvature of 
the Moon. Micrometeorite release, in contrast, is ex- 
pected to be isotropic over the surfece of the Moon, at 
least within a factor of two. 

The ejected products are assigned an initial velocity 
from the surface drawn from the distribution function 
appropriate to the release mechanism. When the parti- 
cle comes back into contact with the surface, there are 
a variety of processes that can occur, introducing an 
array of interesting physics questions. These are inves- 
tigated in the simulations presented here. When the 
atmospheric particle reencounters the planet, it may 
stick to the surface. It may adsorb to the surface long 
enough to partially or fully thermaUze to the local sur- 
face temperature and then be re-emitted. Or it may 
rebound on contact retaining all or most of its incident 
energy. 

The energy exchange at the surface for particles 
that return to the surface and are re-emitted is param- 
eterized by a thermalization coefficient (w) and a con- 
servation coefficient (f) that governs the energy ex- 
change between the particle and the surface. The 
weights applied to v„ a velocity from the Maxwellian 
distribution at the local surface temperature (Aermal 
accommodation) and to Vj the incident particle velocity 
(rebound) total unity. The inbound and thermal veloci- 
ties are added in quadrature with appropriate weights 
to compute the outbound velocity. The conservation 
coefficient is applied afterward to provide a separate 
means of damping particles. If a particle is re-emitted, 
the direction of release occurs with an isotropic angu- 
lar distribution. However, given the microstructure of 
the regolifri, this is a simplification. The re-emitted 
particle is followed on all of its ballistic hops until it is 
lost from the system either to sticking, escape, or pho- 
toionization. 

When the particle encounters the surface, the code 
determmes whether the particle will stick or be re- 
emitted depending on the sticking functions assigned 
to the simulation. In these simulations, sticking is ap- 
plied as a temperature-based function. For each time 
the particle comes into contact with the surface, the 



local surface temperature is queried. If the temperature 
is below the setting, the particle sticks. The code can 
immediately consider the later reemission of a stuck 
particle. For example, when nightside sticking is en- 
abled, the code can assume that the particle is re- 
emitted at the dawn terminator by a specified release 
process (thermal desorption or photon-stimulated 
desorption). This way, one can follow a particle until 
it is lost fi-om the planet by escape or photoionization 
rather than just recycling to the regolith. 

The probability of photoionization or photodisso- 
ciation during a given hop is based on the photoioniza- 
tion time [4-5] and the time of flight in sunlight. The 
particle escapes the simulation when it crosses a prede- 
termined boimdary. Here, we use the Hill sphere (35 
Rmooo) as the boundary for escape. 

Conclosion: Of particular interest are the effects 
of sticking at the surface, thermal accommodation, 
temperature dependence of the yields of photon- 
stimulated desorption, and potential priming of the 
surface by magnetospheric plasma bitting the moon 
during plasma sheet crossings. Data fium the Limar 
Prospector spacecraft, which was operating at the time 
of these observations, will be used to determine the 
times of plasma sheet crossings by the moon. Tem- 
perature-based sticking in the model creates a dawn 
enhancement in the distribution of sodium in the at- 
mosphere. However, compared to data, the density is 
too high if all of the particles that stick to the cold sur- 
faces are released as the surface rewarms. Therefore 
we introduced a reemission percentage to the model. 
Here, we assume a portion of the sodium diffuses back 
into the regolith when the temperature warms, while 
the rest diffuses out of the regolith. This decreases the 
density of the sodium exosphere. 
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